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^ ■ Abstract 

H ; 

5^ I We analyze the forward-backward asymmetry of the decays B (vr, K)i^i with 

^ = /X, or T in the framework of the constrained minimal supersymmetric standard 
model. We find that, the asymmetry is enhanced at large tan and depends strongly on 
the sign of the fi parameter. For /i > 0, the asymmetry is typically large and observable 
whereas for /i < 0, it changes the sign and is suppressed by an order of magnitude. 
Including cosmological constraints we find that the asymmetry has a maximal value 
of about 30 %, produced when Higgs- and gauge- induced flavor violations are of 
comparable size, at a value of tan/3 ~ 35. The present constraints from the B- 
factories are too weak to constrain parameter space, and the regions excluded by them 
are already disfavoured by at least one of BR(i? Xs^), g — 2, and/or cosmology. 
The size of the asymmetry is mainly determined by the flavor of the final state lepton 
rather than the flavor of the pseudoscalar. 



1 Introduction 



There are sound theoretical and experimental reasons for studying flavor-changing neutral 
current (FCNC) processes. Such transitions, being forbidden at tree level, provide stringent 
tests of the standard model (SM) at the loop level. Besides, FCNCs form a natural arena 
for discovering indirect effects of possible TeV-scale extensions of the SM such as supersym- 
metry. Among all the FCNC phenomena, the rare decays of the B mesons are particularly 
important as many of the nonperturbative effects are small and under control. 

In addition to having already determined the branching ratio of -B — > Xg'j [|l|] and the 
CP asymmetry of -B ^ J/tpK [Q, experimental activity in B physics, has begun to probe 
FCNC phenomena in semileptonic B decays p, |^, ^, Therefore, with increasing data 
and statistics, these experiments are expected to give precise measurements on long- and 
short-distance effects in semileptonic decays, e.g. B — > P£~^i~ [P = K, vr) and B — ^ Vi~^i~ 
(y = K*,p). The key physical quantities that can be measured are the branching ratios, 
CP asymmetries and several lepton asymmetries. 

In searching for the physics beyond the SM it is often necessary to deal with quantities 
that differ significantly from their SM counterparts. This is because there are large uncertain- 
ties coming from the hadronic form factors making it hard to disentangle new physics effects 
from those of the hadronic dynamics. For this reason, the pseudoscalar channel B Pi~^i~ 
provides a unique opportunity, since the forward-backward asymmetry Apg in this channel 
is extremely small in the SM (due to a suppression of order m^mfe/M^), and this remains 
true in any of its extensions unless a scalar-scalar type four-fermion operator such as 

= ^VttV*Cimt,)qEbR-U (1) 

provides a significant contribution to the decay amplitude. Clearly, such operator struc- 
tures can arise only from the exchange of a scalar between the quark and lepton lines with 
flavor-violating couplings to the quarks. For instance, by extending the SM Higgs sector 
to two SU(2) doublets, operator structures of the form (|l]) can be generated [0 excluding 
the possibility of ad hoc tree level FCNCs. Although the coefficient C(mb) in eq.(|l|) is still 
proportional to the lepton mass, it can receive an enhancement when the ratio of the two 
Higgs vacuum expectation values, tan /5 is large. 

Supersymmetry (SUSY), is one of the most favoured extensions of the SM which stabilizes 
the scalar sector against ultraviolet divergences, and naturally avoids the dangerous tree 
level FCNC couplings by coupling the Higgs doublet Hu (Hd) to up-type quarks (down-type 
quarks and charged leptons). The soft-breaking of SUSY at the weak scale generates (i) a 
variety of new sources for tree level flavor violation depending on the structure of the soft 
terms, and (ii) radiatively generates various FCNC couplings even if the flavor violation is 
restricted to the CKM matrix. The first effect, which cannot be determined theoretically, is 
strongly constrained by the FCNC data [§], and therefore, as a predictive case, it is convenient 
to restrict all flavor-violating transitions to the charged-current interactions where they 
proceed via the known CKM angles. This is indeed the case in various SUSY-breaking 
schemes where hidden sector breaking is transmitted to the observable sector via flavor-blind 
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interactions e.g. gauge-mediated and minimal gravity-mediated scenarios. This minimal 
flavor violation scheme adopted here is well motivated by minimal supergravity in which all 
scalars receive a common soft mass, mo, at the unification scale. 

The common origin for scalar masses is one of the parameter restrictions which define 
the constrained version of the supersymmetric standard model (CMSSM). The low energy 
sparticle spectrum in the CMSSM is specified entirely by four parameters, and one sign. In 
addition to tuq and tanjS, the remaining mass parameters are the gaugino masses and super- 
symmetry breaking trilinear mass terms. These too, are assumed to have common values, 
mi/2 and Aq, at the unification scale. In principle, there are two additional parameters, the 
Higgs mixing mass, /x, and the supersymmetry breaking bilinear mass term, B, but since it 
is common to choose tan /3 as a free parameter and since we fix the sum of the squares of the 
two Higgs vevs with Mz, these two parameters are fixed by the requirements of low energy 
electroweak symmetry breaking. One is left simply with a sign ambiguity for /i. Therefore, 
the parameters which define a CMSSM model are {mi/2,mQ, AQjtan (3 and sgn{fi)}. 

Although flavor violation is restricted to the CKM matrix, radiative effects still generate 
FCNC transitions among which those that are enhanced at large values of tan (3 are partic- 
ularly important as the LEP era ended with a clear preference to large values of tan/5 0. 
Indeed, it is known that there are large tan /5-enhanced threshold corrections to the CKM 
entries [1^ allowing for Higgs-mediated FCNC transitions For instance, the holo- 



morphic mass term for down-type quarks HdQD'^ acquires a non-holomorphic correction 
HIQD'^ where the latter term proportional to tan/?/167r^, which is not necessarily small at 
large tan /5. 

In what follows, we will compute the forward-backward asymmetry oi B (vr, K)i~^i~ 
decays in the MSSM. After deriving the scalar exchange amplitudes (P, we discuss several 
theoretical and experimental issues and then identify the regions of SUSY parameter space 
for which the asymmetry is enhanced. We compare our results with existing experimental 
and cosmological constraints. 



2 B ^ (tt, K)l^l in Supersymmetry 



In general, the semileptonic decays B — > (vr, K)l^l proceed via the quark transitions h 
{s,d)i~^i~. The decay amplitude has the form 



A 



TT 



(2) 



where ky = —{2mi,/ q^)qv with q^ = sM^ being the dilepton invariant mass. The Wilson 
coefficients, C-j.Cf^, and Cio have been computed to leading order in |]12|. Higher order 0{as) 
corrections, which are available for small s in the SM will not be considered. The 
coefficients C and C will be discussed below. 
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The kinematical range for the normahzed dilepton invariant mass in terms of the lepton 
and pseudo scalar masses is 4mj/M^ < s < (1 — Mp/MbY which includes the vector char- 
monium resonances J/ip^ ip' ^ ■ ■ ■ whose effects are included in the Cg^^{mb, s). Moreover, 
the four-fermion operators for the light quarks develop nonvanishing matrix elements, and 
these are also included in Cg^^ {nib, s). At higher orders in a^, these effects contribute to 



C^^^ {nib) as well ||T3 |. 



The electromagnetic dipole coefficient Cj^^{nib) is contributed by graphs with the W 
boson, charged Higgs, and chargino penguins. The chargino contribution increases linearly 



with tan (3 at leading order [|12| , and the inclusion of SUSY threshold corrections strengthens 
this dependence IH^. This coefficient is directly constrained by the B — > decay rate, and 
the experimental bounds can be satisfied with a relatively light charged Higgs at very large 
values of tan/3. On the other hand, the coefficient of the vector-vector operator Cg^^ {nib, s) 
is generated by box diagrams, and carries a long-distance piece coming from the matrix 
elements of the light quark operators as well as the intermediate charmonium states 



Finally, the coefficient of the vector-pseudovector operator Cio(m;,) is generated by box 
graphs and is scale independent. Both coefficients C^^^ {mb, s) and CiQ{nib) are less sensitive 
to tan/3 than is C^^^{nib). 

Within the SM, these coefficients typically have the values C^^^ {nib) ~ —0.3, C^^^ {nib, s) ~ 
4.4 (excluding its long-distance part), and Cio{mb) ~ —4.7 [|15[ which, however, are allowed 
to vary considerably within the existing bounds The inclusion of SUSY contributions, 
for instance, implies large variations in Cj^^ {nib) (even changing its sign), and typically a 
~ 10% variation in Cl^^{nib,s) and Cio(mfc) ||T7| . 

The scalar-scalar operators in the decay amplitude are generically induced by the ex- 
change of the Higgs scalars and suffer invariably from the nrigmb/M^ suppression. Therefore, 
these operators are completely negligible in the SM. However, in the MSSM, this suppres- 
sion is overcome by large tan (3 effects where the charged Higgs-top diagram is propor- 
tional to tan^/5, and the chargino-stop diagram is ~ (9(tan^/?). In more explicit terms, 
C{mb) = -Cio{mb), and 



2mbmeGF 1 1 

C{nib) = 



v/2 47ra (1 + tan/3)(l + (e^ + /i^e^) tan/3) 
tan^ l3f{xtH) - tan^ ^^^"^ Ap-M^ {f(^x^h) - fi^hh)) 



(3) 



where is the sign of the /x parameter, x is the lighter chargino, Xij = m^/mj, f{x) = 
xloga;/(l — x), and the parameters and eh, which are typically O{10~^), are defined in 
I TTI , [T^ . Finally, ht is the top quark Yukawa coupling, Mj^ is the light stop mass, and At is 
the low energy value of the SUSY breaking top- Yukawa trilinear mass term obtained from Aq 
by the running of the RGBs. Clearly, the charged Higgs contribution, which is the dominant 
one in two-doublet models 0, is subleading compared to the chargino contribution. The 
sign of C{nib) depends explicitly on e^. Therefore, the forward-backward asymmetry in 
B — > (vr, K)i^i'' decays depends strongly on the sign of the /i parameter. 

From the experimental point of view, it is useful to analyze the normalized forward- 
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backward asymmetry defined as 



J~l "-^dzds + JO ^^dzds 



where z = cos 6, 9 being the angle between the momenta of P and 
gives the exphcit expression 



(4) 



'. A direct calculation 



(5) 



where 



+ ste[iB^ois)-Ais)y + v{syAis) 
+ 2t,{l - tp/A - s)[B,^{s) - A{s)]A,^{s 



+ t£tpAio(s) 



(6) 



Here = Amj/M^, tp = AM^/M^, X{s) = (1 - s - tp/Af - stp, v{s) = (1 - tf/s)^/^ ^nd 



Alois) 
A{s) 



Cl^f{m,,s)U{s)-Cf^{m,)Ms), 
Cio{mt)f+{s) , fiio(s) = Cio(m;,)(/+(s) + 



C{mb) 



^ 1(1- Ml/Ml)/. W + .^/-Wl 



(7) 



where fyis) = {2mb) / {Mp + Mp)/r(s). The form factors /+, /_ and fx are not measured at 
present and one has to rely on theoretical predictions. In what follows we use the results of 
the calculation of these form factors from QCD sum rules for both B ^ K and B ^ n 
transitions. 

In general, the hadronic form factors are uncertain by ~ 15%, and this translates into 
an uncertainty of approximately 35% in the branching ratio. Especially for low s, below 
the charmonium resonances, the theoretical prediction for the branching ratio contains large 



uncertainties [0. Therefore, theoretically the large dilepton mass region is more tractable. 
On top of the form factor uncertainties, there are further problems in treating the contribu- 
tions of the charmonium resonances (embedded in the Wilson coefficient C^^^ {nib^ s)). For 
instance, the recent BELLE experiments 0, subtract such resonance contributions by veto- 
ing the range 0.322 ^ s ^ 0.362. Then the experimental bound on the branching ratio turns 
out to be 

0.0328 < 10^ X BR(5 ^ Kfx+fx-) < 2.395 at 90% C.L. (8) 

which we will take into account in making the numerical estimates below. It can be also 
noted that the same decay mode has not been observed by BABAR: BR(i? — > Kfi~^fi~) < 
4.5 X 10~^ at 90% C.L. In addition, the vector kaon final states have not been observed 
yet: BR{B K*ij,+i^-) < 3.6 x lO'^ at 90% C.L. [|, |. One notes that, the asymmetry is 
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s 



Figure 1: The dependence of Apj^{KT^r ) on the normalized dilepton invariant mass, 
s = /M'^, for various (allowed) points in the SUSY parameter space. Each curve is labeled 
as tan {3 where (mi/2, mo) = (400, 100), (400, 100), (900, 700), (1200, 600) GeV from top to 
bottom. One notices that the asymmetry is typically small for fi < as was noted at various 
instances before. 



large in regions of the parameter space where the branching ratio is depleted, and therefore, 
the BELLE lower bound on i? — > Kfi^fi^ is an important constraint which can prohibit the 
asymmetry taking large marginal values. Clearly, in the presence of C{mh), which can take 
large values in SUSY, the would-be experimental constraints on the CiQ{mij)-Cg-^-^ {rrih, s) 
plane are lifted. 

Furthermore, the pure leptonic decay modes, Bg^d depend directly on the Wilson 

coefficients Cio{mb), C{mb) and C{mh). In the SM, BR(Bs — >• ^'^^~) ~ 10~^ which is 
approximately three orders of magnitude below the present bounds BR{Bs /U^/W ) < 
2.6 X 10~^ [0. The SUSY contributions, especially at large tan/?, can enhance the SM 
prediction typically by an order of magnitude, and the bounds can even be violated in certain 
corners of the parameter space |21| . In what follows the constraints from Bs as well 



as the muon g — 2 (as they are directly correlated [^) will be taken into account. We will 



refer to the constraints from B {K , K*) n~ and B^^d collectively as 5-factory 

constraints. 

For the constraints on the SUSY parameter space to make sense it is necessary to be far 
from the regions of large hadronic uncertainties, and thus, below we will restrict the range 
of s to lie well above the charmonium resonances and well below the kinematical end point. 
In Fig. |I], we show the variation of the asymmetry with the normalized dilepton invariant 
mass, s = /M'^ for various values of the SUSY parameters (see below for further discussion 
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of these choices). The irregularities in the s-dependence of the asymmetry are similar to 
those in the B K*i~^i~ decay. The various bumps and valleys come from the relative sizes 
of individual terms contributing to Apg. It should be noted that in the region around the 
value s = 0.75, the s dependence of the asymmetry is rather smooth. Therefore, in forming 
the constant asymmetry contours in the space of SUSY parameters we will take s = 0.75 
(corresponding to = 20.1 GeV^). 

At low values of the asymmetry one should also take into account the final state elec- 
tromagnetic interactions. Indeed, photon exchange between the lepton and P = K,tt lines 
is expected to induce an asymmetry O^a/n) implying that only asymmetries Apg larger 
than ~ 1% can be trusted to follow from SUSY effects, unless the interplay with the electro- 
magnetic corrections is explicitly taken into account. For large values of asymmetry, when 
the observation of the effect becomes feasible, higher order QCD effects (not yet calculated) 
can in principle modify our results somewhat. However, it is highly unlikely that these 
corrections will dramatically reduce the asymmetry discussed here. 

One should also note that in the limit of exact SU{3) flavor symmetry, the asymmetries 
in B ^ Ki^i^ and B ni^i^ decays must be the same. Due to SU{3) breaking effects, 
which show up in different parameterizations of the form factors /+, /_ and fi for B ^ n 
and B ^ K transitions |]18|, their asymmetries are expected to differ slightly. Clearly it is 
the lepton flavor that largely determines the size of the asymmetry, rather than the flavor of 
the final state pseudoscalar. 

Finally, before starting to scan the SUSY parameter space, it is worth while discussing 
the sensitivity of Apg to some of the parameters. First of all, as tan (3 grows, two of the 
Wilson coefficients, Cj-^^ {nib) and C{mf,), grow rapidly up to the bounds obtained from rates 
of the decays {B — >• Xg'j) and {B Kfi~^fi~). Since Re[Cg-^^ {nib, s)] > and Cj^^ {nib) > 0, 
Re[y479(s)] increases with tan/5. However, this increase is much milder than the tan^/5 
dependence of C{nib) causing A{s) to take large negative (positive) values for yU > (/i < 0). 
Therefore, large tan/3 effects influence not only the numerator of Eq. (^ but also the 
denominator S(s) (proportional to the differential branching fraction) via the destructive 
(constructive) interference with -Bio(s) {Cio{mb) remains negative in SUSY) for /i > (/z < 
0). However, as tan/3 keeps growing, depending on the rest of the SUSY parameters, the 
effect of C{nib) eventually becomes more important, and the asymmetry falls rapidly due 
to the enhanced branching ratio. In this sense, the regions of enhanced asymmetry depend 
crucially on the sign of the fi parameter and the specific value of tan /?. Moreover, as the 
expression of C{nib) makes clear, there can be sign changes in the asymmetry in certain 
regions of the parameter space due to the relative sizes of the masses of the lighter chargino 
and stops. Such effects will also give small asymmetries just like the fi < case. 

Our work extends a previous analysis of this asymmetry ||23[ by including the large gluino 
exchange effects (contained in the quantity in Eq. (H)) and the explicit dependence on 
the sign of the /i parameter. In addition, we go beyond the work in as well as in the 
preceding work by resumming the higher order tan /? terms which increases the validity 
of the analysis at large values of tan /? . We note that a computation of the large tan /3 



effects can be carried out in the gaugeless limit |jTT| which eliminates some of the diagrams 
considered in 
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In the numerical analysis below, we will analyze the forward-backward asymmetry of the 
decays B Pi~^i~ by taking into account the above-mentioned constraints from B factories 
as well as other collider and cosmological constraints. We will be searching for those regions 
of the SUSY parameter space in which the asymmetry is enhanced. In particular, we will 
be particularly interested in the sensitivity of the asymmetry to tan/?, the sign of the n 
parameter, as well as the common scalar mass mo and the gaugino mass mi/2- 



3 Results 

In our analysis, we include several accelerator as well as cosmological constraints. From 



the chargino searches at LEP ||2^, we apply the kinematical limit m^± > 104 GeV. A more 
careful consideration of the constraint would lead to an unobservable difference in the figures 
shown below. This constraint can be translated into a lower bound on the gaugino mass 
parameter, mi/2 and is nearly independent of other SUSY parameters. The LEP chargino 
limit is generally overshadowed (in the CMSSM) by the important constraint provided by 
the LEP lower limit on the Higgs mass: rrih > 114.1 GeV f^. This holds in the Standard 
Model, for the lightest Higgs boson h in the general MSSM for tan/5 <^ 8, and almost always 
in the CMSSM for all tan/3. The Higgs limit also imposes important constraints on the 
CMSSM parameters, principally mi/2, though in this case, there is a strong dependence on 



tan p. The Higgs masses are calculated here using FeynHiggs p6[, which is estimated to 
have a residual uncertainty of a couple of GeV in m/j. 

We also include the constraint imposed by measurements of 6 57 [jl], These agree 
with the Standard Model, and therefore provide bounds on MSSM particles, such as the 
chargino and charged Higgs masses, in particular. Typically, the 6 — 57 constraint is more 
important for < 0, but it is also relevant for /i > 0, particularly when tan/3 is large. 

The final experimental constraint we consider is that due to the measurement of the 
anomalous magnetic moment of the muon. The BNL E821 [27] experiment reported a 
new measurement of = \{g^ — 2) which deviates by 1.6 standard deviations from the 
best Standard Model prediction (once the pseudoscalar-meson pole part of the light-by-light 
scattering contribution pSf is corrected). Although negative values of n are no longer entirely 



excluded [^, the 2-a limit still excludes much of the /i < parameter space [^. < is 
allowed so long as either (or both) mi/2 and mo are large. 

We also apply the cosmological limit on the relic density of the lightest supersymmetric 
particle (LSP), = VLy^PcriUcaU and require that 

0.1 < n^h^ < 0.3 (9) 

The upper limit is rigorous, and assumes only that the age of the Universe exceeds 12 Gyr. 
It is also consistent with the total matter density Vl^ ^ 0.4, and the Hubble expansion rate 
h ~ 0.7 to within about 10 % (in units of 100 km/s/Mpc). On the other hand, the lower 
limit in (|^) is optional, since there could be other important contributions to the overall 
matter density. 



The cosmologically allowed regions in the CMSSM have been well studied [p^, p2|. There 



are generally large, 'bulk' regions of parameter space at low to moderate values of mi/2 and 
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mo at all values of tan (3. There are additional regions which span out to large values of 
mi/2 due to co-annihilations with light sleptons, particularly the lighter f |Q. At large 
tan /3, there are also regions in which the lightest neutralino sits on the s-channel pole of 



the pseudo-scalar and heavy scalar Higgs producing 'funnel'-like regions [34, 31]. Finally 



there are the so-called 'focus-point' regions ||3^ which are present at very large values of mo- 
Generally, these regions have a lower asymmetry (because of the large value of mo), however, 
at values of tan/3 ~ 50, asymmetries as large as 10% are possible. 

In Fig. we show the contours of constant Apg(i^'r+r~) in mi/2^tan/3 plane for 
Aq = 0, mo = 100 GeV and /i > 0. The constraints discussed above are shown by various 
curves and shaded regions. The nearly vertical dashed (black) line at the left of the figure 
shows the chargino mass constraint. Allowed regions are to the right of this line. The dot- 
dashed Higgs mass contour (red) labeled 114 GeV, always provides a stronger constraint. 
Allowed regions are again to the right of this curve. However, one should be aware that there 
is a theoretical uncertainty in the Higgs mass calculation, making this limit somewhat fuzzy. 
The light (violet) solid curve shows the position of the 2-a g — 2 constraint, which again 
excludes small values of mi/2- In the dark (red) shaded region covering much of the upper left 
half of the plane, the lighter f is either the LSP or is tachyonic. Since there are very strong 
constraints forbidding charged dark matter, this region is excluded. The medium shaded 
(green) region shows the exclusion area provided by the 6 — 57 measurements. Finally the 
light (turquoise) shaded region shows the area preferred by cosmology. Outside this shaded 
region, the relic density is too small and is technically not excluded. 

Putting all of the constraints together, we find that for this value of mo = 100 GeV 
and /i > 0, the allowed region is bounded by 300 GeV ^ mi/2 ^ 500 GeV and 5 ^ 
tan /3 ^ 20. In the allowed region, the forward-backward asymmetry varies rapidly from 
very small (unobservable) values up to 10%. There is a wide region with an observable 
1-10% asymmetry though the 10% region is quite narrow (restricted to tan/? ~ 20). 

In Fig. 03, we show the corresponding result for the opposite sign of fi. While the 
cosmologically allowed region is qualitatively similar to the fi > case and the Higgs limit 
is slightly stronger, we see that the b ^ sj constraint is significantly stronger. Indeed, the 
combined constraints from b ^ sj and a f LSP exclude tan /5 > 15 for this value of mo and 
fi < 0. The 2-0" constraint from (7 — 2 is also significantly stronger and when combined with 
the f LSP constraint now exclude values of tan /3 > 8. 

As mentioned earlier, for /x < both the sign and size of the asymmetry have changed. 
In general, the size of the asymmetry is suppressed by an order of magnitude. Clearly, the 

6 — >■ 57 constraint now allows only a small region with a —0.3 1% asymmetry. However, 

when all constraints are combined they exclude almost completely the otherwise allowed 
regions. At higher values of mo, slightly larger asymmetries are possible. At mo = 200 GeV 
(with /i < 0), 6 ^ 57 allows asymmetries as large as —2%, however, the g — 2 data still 
restricts the asymmetry to values below about —0.4 %. Even at large tan/? and very large 
mo, we will see below that for < 0, asymmetries never exceed ~ —1 %. We note that 
independent of the sign of fi, the asymmetry is maximized for intermediate values of tan (3, 
i.e. it does not monotonically increase with increasing tan/3 as was already argued earlier. 
The main conclusion from this figure is that the sign of the fi parameter must be positive in 
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Figure 2: The constant A-p'q{Kt^t ) contours in mi/2^tan/9 plane for Aq = 0, with 
mo = 100 GeV, /i > (a); ttiq = 100 GeV, /i < (b); mo = 200 GeV, /i > (c); and 
mo = 300 GeV, fx> (d). We take mt = 175 GeV and mb(mb)|^ = 4.25 GeV. In all of 
the panels the black dashed line shows the 104 GeV chargino mass contour, the dot-dashed 
(red) curve stands for mt = 114 GeV (evaluated using the FeynHiggs code f^). The light 
(turquoise) areas are the cosmologically preferred regions with 0.1 < Vl^h'^ < 0.3. In the 
dark (brick-red) shaded areas the LSP is Ti and thus excluded. The medium (green) shaded 
regions are excluded by BR{B — >■ ^^7). The light (violet) curve shows the — 2 (the area 
to the right of which is allowed). 
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order to have large observable Apg(ii'r+r^). 

We note that there are already 5-factory constraints due to recent BELLE and BABAR 
experiments [0, ^ |[. For > 0, they exclude a small region (not plotted) with mi/2 ^ 
120 GeV and tan/5 > 43.5 (lying in the region with a charged LSP), whereas for n < the 
excluded region is shifted to mi/2 ^ 260 GeV and tan/3 > 30.5 (now lying in the region also 
excluded by BR(5 ^ X,7)). 

The behaviour observed in Fig. ^) and b) is subject to large variations once the GUT- 
scale input parameters are varied. This is seen in Fig. ^) and d) where the constant 
asymmetry curves are plotted in mi/2-tan /3 for ^ > 0, Aq = 0, and tuq = 200 GeV and mo = 
300 GeV respectively. It is clear that with increasing mo the cosmologically preferred 'bulk' 
region is shifted towards larger tan/3 values making it possible to get larger asymmetries. 



In addition, we see very clearly the effect of x "~ coannihilations |Q which extend the 
cosmological region to high values of mi/2- The region below the 'bulk' and coannihilation 
region is excluded as it corresponds to an area with Qh"^ > 0.3. While the chargino, Higgs, 
and g — 2 constraints are only slightly altered at the higher value of mg, we see that the 
charged LSP constraint is relaxed in Fig. |^) and greatly relaxed in panel d). 

For the higher values of rriQ we see from Fig. |]c) that for 300 GeV <^ ^ 800 GeV 
and 25 ^ tan/? ^ 40 the forward-backward asymmetry ranges from 1% to 30% in regions 
which are not excluded by any experimental or cosmological constraints. In particular, when 
350 GeV ^ mi/2 ^ 400 GeV and 30 ^ tan/3 ^ 35 the asymmetry is well observable with a 
typical ~ 30% peak value. Here the S-factory constraints are effective for mi/2 ~ 100 GeV 
and tan (3 > 47.5, i.e. only a small region in the upper left corner. 

For the tuq = 300 GeV case shown in Fig. ^), the cosmologically allowed region is 
now shifted up past the maximum of Apg(i^rr), and is now typically 10%. Overall the 
forward-backward asymmetry is larger than 1% in the cosmologically allowed region and 
extends over the range 300 GeV ^ 1711/2 ^ 1000 GeV and 25 ^ tan/3 ^ 50. 

We have also checked some cases with nonzero values of Ao, assuming it to be either 
constant {e.g. set to 2 TeV), or varying in proportion with mi/2 {e.g., Aq = 21711/2). For 
a variable Aq, results were not qualitatively different from the results shown here. For a 
large and fixed value of Aq, the comso logical regions of interest could be very different [^|, 
however, the asymmetry was found to be quantitatively similar as the results shown here. 
However, we can not claim to have made a systematic examination of the Aq ^ parameter 
space. 

In Fig. II, we show the contours of Apg(7rr+r^) (left panel) and Apg(i^/i"'"/i~) (right 
panel) for Aq = 0, tuq = 200 GeV and /i > 0. A comparison of the left panel with panel 
c) of Fig. ^ shows that there is very little difference between 7ri~^i~ and K£~^i~ final states 
as far as the asymmetry is concerned. Indeed, as mentioned before, the difference between 
the asymmetries is a measure of the SU(3) flavor breaking or the different parameterizations 
of the associated form factors. Therefore, the similarity or dissimilarity of these two figures 
depends on how the hadronic effects are treated for the kaon and pion final states. On the 
other hand, the comparison between panel c) of Fig. ^ and the right panel of Fig. |^ shows 
that the asymmetry is suppressed for the f7'^fi~ final states. The asymmetry does not reach 
the 1% level in any corner of the allowed regions. 
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mo = 200GeV, |i > niQ = 200 GeV, |i > 




' I I I I""' I I I I I ""1 = I I I \'" I I I I I ""1 
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ini/2 (GeV) mi/2 (GeV) 



Figure 3: The contours of A-pj^i^nr^r ) (left panel) and A-pj^^K fi'^ fi ) (right panel) for 
Ao = 0, mo = 200 GeV and /x > 0. All other curves and shaded regions are taken from Fig. 

In Fig. we show the Apg(_ft'r"'"r~) contours as projected onto the mi/2-mo plane 
for Aq = 0, tan (3 = 10, /i > 0. The curves and shading in this figure are as in Fig. ^ Seen 
clearly are the 'bulk' and coannihilation cosmological regions. The latter traces the border 
of the f LSP region which is now found at the lower right of the figure. The constraints from 
b ^ s'y and g — 2 only exclude a small portion of the parameter plane at the lower left. In 
this figure, it is the area above the cosmological shaded region which is excluded due to an 
excessive relic density. 

This figure illustrates the dependence of the asymmetries on the chargino and stop masses, 
and indeed, the asymmetry changes sign at fairly large mo and small mi/2. This effect results 
in a sign change in C{mh) due to the competition between the stop and chargino masses. As 
Fig. ^) makes clear, the asymmetry is positive and remains below the 0.3% level yielding 
essentially no observable signal at all. 

An example with tan/5 = 10 and < is shown in Fig. ^). While the cosmological 
region is similar to that in panel a), the constraint from 6 — > 37 is significantly stronger, as 
is the constraint from g — 2. The latter imply that the asymmetry is never much larger than 
—0.1%. As one can see, it is difficult to get any observable effect at low values of tan /5, for 
any sign of the fi parameter. 

In the lower two panels of Fig. ^, we show the contours of Apg(i^r^r^) for larger values 
of tan/3. As one can see, the 'bulk' cosmological regions are pushed to higher values of 
mo and we also see the appearance of the 'funnel' regions where the LSP relic density is 
primarily controlled by H,A s-channel annihilations. In panel c), there are large positive 
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tan P = 10 , )i > tan P = 10 , |i < 




mi/2 (GeV) mi/2 (GeV) 



Figure 4: The constant Apj^{KT^T^) contours in mi/2^rnQ for Aq = 0, and tan (3 = 10, 
/i > (a); tan/? = 10, fi < (h); taii/S = 50, /z > (c); and tan/3 = 35, < (d). 
Contours and shaded regions are as in Fig. ^ 
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asymmetries (> 10%) in a broad region extending from (mi/2,mo) = (300,300) GeV all 
the way up to (7Tii/2, ""^o) = (l-S? 1-1) TeV. The 5-factory constraints only exclude a small 
region in the bottom left corner bounded by mo ^ 200 GeV and mi/2 ^ 120 GeV. Panel 
d) of Fig. ^, on the other hand, shows that, for negative /i and large tan/5, the allowed 
range of asymmetry is around —1% never reaching the —2% level mainly due to the — 2 
constraint. There are three cosmologically preferred strips, the wider one located in the 
region (mi/2,?^o) ^ (900,700) GeV where the asymmetry is at most —1%. For /i < the 
5-factory constraints exclude two small regions bounded by uiq <^ 200 GeV and mi/2 ^ 
280 GeV, and mo ^ 200 GeV and mi/2 ^ 1940 GeV which both are already excluded by 
BR{B Xs'j) and cosmology. 

4 Summary 

We have discussed the forward-backward asymmetry of B ^ {K, TT)i^i~ decays which is 
generated by the flavor-changing neutral current decays mediated by the Higgs bosons. 
In addition to the known properties e.g. the smallness of the muon production asymmetry 
compared to the r production, the approximate independence of the asymmetry to the flavor 
of the final state meson, the reduction of the hadronic uncertainties in the high dilepton mass 
region, etc. we find that 

• The remarkable enhancement of the asymmetry is a unique implication of SUSY not 
found in the SM and its two-doublet version. 

• The regions of large asymmetry (> 10%) always require /i > 0, and when /i changes 
sign so does the asymmetry with an order of magnitude suppression in its size. 

• The asymmetry in the decay channels with the lepton pair r+r^ is significantly larger 
(approximately by the factor m^/m^) than in the channels with Thus in spite 
of a greater difficulty of restoring the kinematics in the events with r leptons, these 
may still be advantageous for measuring the discussed asymmetry. 

• The asymmetry is not a monotonically increasing function of tan/? instead it is max- 
imized at intermediate values above which the scalar FCNC effects dominate and en- 
hance the branching ratio, and below which such FCNC are too weak to induce an 
observable asymmetry. 

• Though BR(i? —>■ Xs'j) strongly disfavors the negative values of fi, the present bounds 
from the muon g — 2 measurement are much stronger, and it typically renders the 
asymmetry unobservably small. 

• The cosmological constraints are generally very restrictive. When mo ~ 200 GeV, 
■mi/2 ~ 400 GeV, tan/3 ~ 30 and /i > there is a relatively wide region where the 
asymmetry is ~ 30% for r, and typically ~ 0.5% for muon final state. 
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• The 5-factory constraints are generally too weak to distort the regions of observable 
asymmetry, and the regions excluded by them are already disfavoured by one or more 
of BR{B —> Xg'j), fi LSP and (7^ — 2. With increasing statistics it is expected that 
the branching ratios of semileptonic modes will be measured with better accuracy so 
that, for instance, the allowed range in (|^) will be narrowed. In this case, there may 
be useful constraints on regions of enhanced asymmetry. 
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